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PHASE DISTORTION DUE TO AIRFLOW OVER 
A HEMISPHERICAL LASER TURRET 



by 

A. E. Fuhs* and S. E. Fuhs** 

I. INTRODUCTION 

The index of refraction of air is related to mass density. Starting 
at modest subsonic Mach numbers, e.g. 0.3, and extending to transonic 
flight, flow over a blunt laser turret causes significant changes in 
density. The change in density, which in turn changes index of refraction, 
can defocus a laser beam. 

Variable index of refraction within context of aeronautical appli- 
cations is not a new problem; however, there are two aspects which are 
new. First, the laser beam is large relative to a characteristic turret 
dimension. Second, for many airborne applications the laser telescope 
requires a turret which protrudes into the flow. 

Reconnaissance aircraft fly photographic missions. The compressible 
flow over the aircraft acts like a distorted lens and conceivably could 
cause distorted photographic images. The variable index of refraction is 
not a serious problem for reconnaissance aircraft for two reasons. The 
aperture of the camera is small compared to the characteristic scale of 
the distortions; in the case of a laser telescope this in not true. Fur- 
thermore, the camera on a reconnaissance aircraft is not mounted in a 
turret. To aim the camera the pilot aims the aircraft. Camera windows 
are flush with the aircraft surface. 

Cruise missiles with intercontinental ranges may use celestial navi- 
gation which requires precise measurement of angular location of stars. 

* Distinguished Professor of Mechanical Engineering 

** Student, California Institute of Technology 
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A distorted lens created by the ambient flow field causes an error in 

measurement. Knowledge of the external aerodynamics is important. 

A model which is useful for understanding the density field in 

the vicinity of a turret is the flow around a sphere or an ellipsoid. 

The plane containing an incoming fluid velocity vector and a diameter 

of the sphere does not have flow components normal to the plane. This 

is a plane of symmetry and conceptually can be considered a wall. The 

flow over the hemisphere above the plane can be calculated readily. 

For incompressible flow, the solution for the flow over a sphere 

is well known. See for example pages 464-465 of Milne-Thompson,^^^ 

pages 92-93 of Lamb or pages 339-342 of Karamcheti.^ Techniques 

are available to introduce the influence of compressibility or, stated 

another way, the influence of Mach number. One method is the Janzen- 

(4) 

Rayleigh technique discussed briefly by Van Dyke. Another brief 
description is given in Chapter 10 of Liepmann and Puckett . ^ ^ ^ The 
Janzen-Raylei^ method is discussed also on pages 328-334 of Oswatitsch. 

(7) 

Lord Rayleigh provides the second order solution for a sphere in a 
paper published in 1916. 
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II. DISCUSSION OF ANALYTICAL MODELS 



Optical distortion due to the external flow field can be divided 
into two categories. Viscous flow phenomena fall into one category and 
include shear layers, laminar and turbulent boundary layers, and shedding 
of discrete vortices. The other category involves the inviscid flow 
field. This paper treats the latter category. 

Laser turrets can be classified as blunt or small-perturbation 
turrets. If the turret has a surface with a normal vector, n, aligned 
parallel, or nearly so, to the freestream velocity vector, the turret 
is "blunt". If the angle between the normal vectors of the turret sur- 
face and the freestream velocity vector are everywhere large, e.g. 60° 
to 90°, then the turret is "small-perturbation". 

The turret being considered here is a hemisphere, which is obviously 
a blunt turret. At the front stagnation point of a hemispherical turret, 

^ • n = - 1.0. 

For a blunt turret, compressibility effects set in at a smaller 
freestream Mach number, M^. A useful concept is the critical Mach number, 
M^*. At the critical Mach number, somewhere on the body a local Mach 
number is just sonic. When M^ > M^*, the flow becomes inherently nonlinear 
and shock waves appear. 

For a circular cylinder with its axis normal to the flow, the critical 

Mach number is 0.3985. For a sphere, the critical Mach number is 0.5868. 

A hemisphere-cylinder, which is oriented with the cylinder axis parallel 

to the freestream flow direction, has a critical Mach nun±)er less than 
( 8 ) 

0.7. Hsieh gives experimentally determined pressure coefficients for 
a hemisphere-cylinder over the Mach number range 0.7 to 1.0. In refer- 
ence 9 Hsieh extends the Mach number range from 1.05 to 1.82. Using 
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from reference 8, the maximum local Mach number is 1.25; the maximum 
local Mach number occurs at an angle of approximately 60° on the sphere. 
The angle is measured from the stagnation point to the location of 



minimum C^. 



For small perturbation turrets the critical Mach number is larger. 

( 10 ) 



(lalculations were made of a turret on a cylindrical fuselage, 
turret shape was determined by 



The 



R(x,e) = Rq + ^ (1 + cos 1^) (1 + cos f6) 
R(x,e) = R^ for I < |x| < L and 7i/f < |0|< n 



( 1 ) 



where R^ is fuselage radius; e, height of turret; £, length of turret; 
1/f, the fraction of the 2n circumference occupied by the turret; and, L, 
the separation between periodic turrets. Table I summarizes the critical 
Mach number. 



TABLE I. 


Critical 


Mach Numbers 


for Small-Perturbation 


Laser Turrets 


£/R„ 


L/R(, 


1/f 


e/R. M * 

0 “ 


maximum slope 
on turret 



1.005 



5.0 



,333 



0.05 


.88 


4.5° 


.10 


.82 


8.9° 


.15 


.76 


13.2° 


.20 


.72 


17.4° 


.25 


.68 


21.3° 


.30 


.65 


25.1° 


.35 


.62 


28.7° 


.40 


.59 


32.0° 


.45 


— 


35.1° 


.50 


— 


38.0° 


.55 


— 


40.7° 
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There are different analytical models that can be developed for 
laser turret geometries. Flow over blunt turrets can be determined using 
expansion techniques such as the Janzen-Rayleigh method mentioned ear- 
lier. As more terms in the expansion are evaluated, the more accurate 
are the flow field calculations as Mach number increases. Upper bound 
is the critical Mach number. For M > M *, the calculations will be 
grossly in error. 

There are numerical techniques which directly integrate appropriate 
differential equations. One method, which is applicable to axisyrametric 
flow, is the technique of South and Jameson. 

For small-perturbation laser turrets the equations can be linear- 
ized. Compressibility effects in subsonic or supersonic flow are incor- 
porated. Transonic flow is excluded since the transonic equations are 
nonlinear even for small perturbation. 

Figure 1 summarizes the regions in a body slope vs. Mach number 
map. The flow within the region ABCD can be described with sufficient 
accuracy using the linearized, small-perturbation equations. The line 
CD represents the critical Mach number. The line BC is shown at a 
slope of 30°. The assumptions involved in the linearization of the po- 
tential flow equation becane less and less valid as the maximum slope 
increases. The region CDJ, which resembles an inverted triangle, requires 
solution of the nonlinear small-perturbation, transonic, potential equa- 
tion. 

The region EFGH defines the area where Janzen-Rayleigh technique is 
most useful. The line GH, which defines the critical Mach number, is 
shown dashed; the critical Mach number is a function of turret shape. 

Hence GH is intended to qualitatively suggest the upper bound for Mach 
number. 
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MAXIMUM BODY SLOPE 



2 



F 



3 4 



e 




FREESTREAM MACH NUMBER , Mqo 



Figure 1. Freestream Mach Number Regions for Application of 
Various Analytical Models. 



There are vertical lines 1, 2, 3, and 4 within region EFGH. These 
lines represent the upper bound of applicability of the first, second, 
"thii^i and fourth order solutions when Janzen- Rayleigh technique is 
used. The first order solution may be applied within region FllE; the 
second order solution, within region F22E; etc. 

A hemispherical turret considered in this paper is an example of 
a blunt turret falling along the line FG. 
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III. JANZEN-RAYLEIGH EXPANSION TECHNIQUE 



In spherical coordinates , the equations of motion are as follows : 



u 



u 



2 2 
^ + 1 9p - a 3p 

3r r 96 r ~ p 9r ~ p 9r 



2 

^ + X ^ + 1 l£-_^ 

9r r 96 r " p r 96 “ pr 96 



( 2 ) 



(3) 



i [u iP + X iP 1 
p |_ 9r r 96J 



1 



9 , 2 . 

^(ru) 



r sin 6 



3 / 

5e 



sin 6) = 0 (4) 



Equations (2) and (3) are Euler's equations. Equation (4) is the con- 
tinuity equation. For flow over the sphere, there is no dependence on 
the azimuthal angle in polar coordinates. Consequently, equations (2) 

and (4) are two-dimensional equations depending only on r and 6. 

2 2 

Multiplying equation (2) by + u/a and equation (3) by - v/a , 
adding and combining with equation (4) yields 



1 9 f 2 . , 

P ^ u) + 



1 9 / ■ 




2 

9u , uv 9v , 9v , V 9v 
9r r 30 3r r 39 



(5) 



The speed of sound, a, is the local value and changes from point to point 
in the flow. It is necessary to account for variations in a. Using the 
energy equation along a streamtube, one can show that 




a 



oo 




2 



V 



( 6 ) 



Inviscid flow is considered. A potential can be introduced for the 
velocity components u ~ U — (7 



7 



( 8 ) 



V 



U 



r 



86 



where U is the freestream velocity at upstream infinity. Introducing 
equations (7) and (8) into equations (5) and (6) yields equations in 
terms of a potential. 

The Rayleigh-Janzen expansion considers the potential function to 
be given as 

<}) = <t>r, + (|), + + (9) 

U‘ CX3^J_ Z ••• 

The solution for (f)^ is considered the first order solution. The solution 
involving both (|)q and (})^ is called the second order solution. In this 
paper only the second order solution is considered. For ease of writing, 
the following definition is introduced: 

ij; = <(), (10) 

To digress briefly from the development of the Janzen-Rayleigh 

expansion, the solution for the flow over a cylinder normal to the fluid 

(4) 

velocity will be discussed. Van Dyke indicates the solution for the 

16 

cylinder has been performed to . For y = 7/5, the result for maximum 
surface speed is 



Snax 

U 



2.00000 + 1.16667M^ + 2.57833m‘^ + 
00 00 

+ 25.59041M^ +....+ 7450M^® 



( 11 ) 



The critical Mach number for a cylinder is = 0.39852 ^ 0.00020. 

In terms of the potential function, the speed of sound becomes 
2 

a 

00 

where F is given by 
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(12b) 



[■- W'?)] 



With substitution of equations (9) and (10), the potential function 
becomes 



v2({, + v^\p = M- 

0 » 



O 11 

(b "b 9 d) ^ U 4* 4’ n n 

^rr r'^ ^0 ^0 ^00 



(13) 



where all terms on the right-hand side are determined from the first 

order potential function (fg. The result of equating coefficients of like 
2 

powers in M yields two equations 



and 



v2<j) = m2 

00 2 oo 



V2(j). = 0 



2 11 
^r*^rr r^ '^r'^0*^0r r^ *^0 '^r'*' r^ *^0 *^00 



(14) 



(15) 



The right-hand side of equation (15) is known due to the fact that 
equation (14) can be solved to yield <J)g. To emphasize this fact, the 
first term on the right-hand side is written explicitly showing ((>g, 

92a 



4^4 
^r rr 






^0 



3r" 



Equation (15) is Poisson's equation. 
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IV. FIRST ORDER SOLUTION 



The first order solution is obtained from Milne-Thompson,^^^ Lamb,^^^ 
(3) 

or Karamcheti. It is 



<}>Q(r, e) = U (r cos 0 + ^ ^ 

where R is the sphere radius. The upstream stagnation point, P, is at 
P(R, 180°); and the downstream stagnation point is at P(R, 0°). Using 
equation (16) one can determine the derivatives to substitute into equa- 
tion (15). These are 



p = U cos 0 (1 - ^ ) 
r r'^ 



rr 



= U cos 0 



= - r U sin 0 



'00 



= - r U cos 0 



'0r 



= - U sin 0 



[‘•iU 

[.•iS] 

[•- 



r3 



The result of combining equations (15) and (17) to (21) is 



(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 



V2<j)^ = |3X(1-X)2 + sin2 0 1^4 - I X + 6X2 _ ^ ^ (22) 

3 3 

where X = R /r has been introduced. 

Before proceeding with second order solution, a comment about 
boundary conditions is appropriate. The component of velocity normal to 
the surface of the sphere must be zero. At infinity, the velocity should 
be U. Consequently, 
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Since 4 >q yields U at infinity, the velocity perturbation due to 
should be zero at infinity. 
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V. RAYLEIGH SOLUTION FOR SECOND ORDER FUNCTION 

(7) 

In July, 1916, Lord Rayleigh reported the solution for the 
second order function. 

The Legendre functions 

= cos 0 (24) 

Pg = j cos^ 0 - y cos 0 (25) 



are introduced into equation (15) in the following manner: First, 

equations (17) through (21) are substituted into equation (15); second, 
terms are grouped so as to form the functions P^ and P^. The result is 





f 36 R6 ^ 


3 \ P + / 




24 r6 ^ 3 


\ p "I 


2[l 


5"F' ^ 


2 \ 


K 

^ ■ 


' TrT^ 2 


rioj ^al 



= 



Equation (26) is in the form 



(26) 



= E E r”"^ P C 



n m 



n nm 



(27) 



where C^ are zero except for m = 7 and 10; all are zero; and C^^^^ 
are zero except for m = 4, 7, and 10. Non- zero C^ are the appropriate 
constants appearing in equation (26). The Legendre differential equation 
is 



^ 

sin 0 d0 



r T 

l^sin e . n 



(n+1) P = 0 
n 



(28) 



Define the function <j> which satisfies the equation 

nm ^ 



r ^ C = v2((, 
n nm ^ 



- Jim ^ ^nm ^ 
nm ” 9r2 r 9r r2sin 



1 9 r • r, ^'^nml 

iSTe 36 « — J 



(29) 



Assume 



nm 



P A 
n nm 



(30) 
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Consequently, 



94> 



„ = (2-m) r PA 

9r n nm 



(31) 



and 



= (2-m) (1-m) r ^ P A 
9r^ n nm 



(32) 



Clombining equations (29) and (32) gives 



C = (2-m) (1-m) r P A 
n nm n nm 



+ 2 (2-m) r'^^P A 



-m+2 



n nm r2 sin 6 



§6 [= 



9P 



sin 6 



n 



96 



nm 



Substitution of equation (28) into equation (33) yields 

r'^PC =A r^P f(2-m) (1-m) + 2 (2-m) - n (n+1) 
n nm nm n I 

As a result 



A = 



nm 



C 

rma 

(2-m) (3-m) - n (n+1) 



(33) 



(34) 



(35) 



Using equations (26), (30), and (35), that part of the second order 



function identified as the particular solution is obtained 






3R3 

10r2 



3R6 

lOrS 



3R5 

17^ 



The homogeneous solution is 



(36) 



(|)^ = Ar P^ + BP^/r2 + Cr^P^ + D P^/r'* 



(37) 



At infinity the potential must be finite and the velocity must be U; 
the function 4 >q of equation (16) yields U. Clonsequently both A and C 
equal zero. At the surface of the sphere 
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=0 



and 



=0 



(38) 




1 

r 



8<t>] 



/9=0 



From the requirements of equation (38) combined with equations (36) 



and (37) 



and 



B = 



3 



U m2 

00 



D = 



27^ m2 



The final expression for the potential function is 






1, 




) 2( 


R3 R6 r5 


r 




?, V + UM 1 


3r2 5r^ 24r8 



+ 



3R3 ^ 27R5 3R6 3R5 

10r2 55r^ " lOr^ 




(39) 



(40) 



(41) 



The first term in brackets is and the second, long term is the second 
order contribution. 
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VI. DEVELOPMENT OF PHASE DISTORTION 



Use of Potential Function to Obtain Index of Refraction 
The index of refraction and density are related by 

n = l + K'^ !1=1 + k— (42) 

Poo PSL P- 

where k' is a constant with value of about 2 x 10~^. Note 
K = flow is assumed to be isentxopic. Consequently 




1 + ^ m2 
2 ” 

1 + ^ m2 



1 

Y-1 



(43) 



The local Mach number is given by 



m2 3 




(44) 



In equation (42) p is the local density, is the freestream density, 
and p^j^ is the sea level density on a standard day. By combining equa- 
tions (6), (12a), (41), (42), (43), and (44), an expression is obtained 



for the index of refraction. 

The turret has a characteristic length; for the sphere this is 
the radius R. The laser beam has a characteristic size which is the 
diameter D or radius R' . The wavelength of the laser radiation. A, 
is also a characteristic length as will be demonstrated in subsequent 
discussion. Two nondlmensional lengths can be formed. Since R occurs 
prominently in the equations for the potential, R is chosen as the 
reference length. With that choice, A/R and D/R are two of the para- 
meters specifying the problem. 
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Phase distortion is a function of the direction that the beam is 
pointed. Define e and a as elevation and azimuth angles. For the 
hemispherical turret, symmetry exists about the diameter of the sphere 
which is parallel to the freestream velocity vector. In general, two 
angles, e and a, are needed to define the beam direction; however, 
for the hemispherical turret only one angle is needed. Looking ahead, 
the angle 6’ occurring in equations (45) and (46) will be used as the 
beam angle; see also Figure 2. 



Geometrical Relations for Calculation of Phase Distortion 

Figure 2 shows the geometry of the problem. The radius of the 
sphere is R. The beam is elevated at an angle 6' . The z axis is the 
same as the polar axis for polar coordinates. 

To locate points within the beam, three variables are used; s is 
distance along the beam, as shown in Figure 3; a is an angle measured 
clockwise from the windward edge of beam; and R' is the radial distance 
from the beam centerline or axis. 

The flow properties, i.e., density and local Mach number, are 
functions of r and 9. These are polar coordinates. Consequently, it is 
necessary to express r and 0 in terms of 9', s, a, and R. The equations 



are 



= j^(R + s) sin 0' - R' cos a cos 0 ' j + R’^ sin 



[‘ 



+ ( R + s ) cos 0 ’ + R ’ cos 



a sin 0’J 



(45) 



and 



9 = tan 






(R + s) sin 0’ - R' cos a cos 9' 



^ + R'^ 



R' sin a 



(46) 



(R + s) cos 0' + R’ cos a sin 0' 



Equations (45) and (46) can be used to evaluate equations (41) to (44). 
The index of refraction is given by equation (42). Define N as the 
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Figure 2. Geometry for Phase Distortion Calculation. 




Figiare 3. Cross Section of the Turret and Beam in the x-z Plane. 
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N A 



number of waves along a distance s. Hence 



dN = P (47) 

A 

where X is the local wavelength. The local wavelength is a function of 
0 ’ , a, R' , and s. It is 



X 




X (1 + k) 



C _ C _ “ 

nf f(l + Kp/p^) ~ (l + Kp/p^) 



As a result, the number of waves is 



(48) 



N = 



X (1 + k) 



/ 



(1 + 



K — )ds 



-S' 



(49) 



A reference is needed. The reference which is convenient is the number 
of waves along the beam axis; designate that quantity as N^. Using this 
reference 



N-N 



1 

0 X (1+k) 



oo 





(1+K— )ds 
P„ 



(50) 



Note the density on the beam centerline for -s'^ s 0 is denoted by p^, 
a reference density. For 0 ^ s < ", the density on the centerline is 
Pq which is a function of s. In equation (50) the value of unity in 
each integrand can be cancelled giving 



,, _ kR 

^"^0 ■ X (1+k) 






P-Pf 




ds 

R 



(51) 



The second integral in equation (51) has been termed the GAP INTEGRAL 
since it enters the equation as a result of the gap between s = -s' and 
s = 0. See Figure 3. In the calculatLors which were conducted, the gap 
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integral was evaluated by assuming p = p ; this is an arbitrary as sump- 
tion. Note that in equation (51) wavelength is nondimens ionalized by R 
to form ^^/R. The distance along a ray is also nondimens ional, i.e., 
ds/R. 

A computer program has been developed for the HP9830 which calculates 
N - Ng as a function of sphere radius, R; elevation angle, 9' ; radius 
within the beam, R' ; angle within the beam, a; freestream Mach number, 

M; wavelength, constant in density equation (42), k'; ratio of heat 
capacities for air, y; freestream speed of sound, a^. Note that 
it is not necessary togjecify p^ since all calculations depend on den- 
sity ratio; see equations (48) and (51). Values used in the calculations 
are as follows: 

R = 0.9144 meter 

"“6 

A =3.8 micron = 3.8 x 10 meter 

00 

k' = 2 X 10"*^ 
a = 342 meter/ sec 
Y = 1.4 

M^, R' , a, and 9' have been varied. Calculations have been performed 
to the extent necessary to plot phase distortion maps. 
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VII. RESULTS OF CALCULATIONS 



Graphical Presentation of Results 

Plots of isocontours of phase shift were made for steps of 18° 
starting with 6' = 0 to 0' = 90°. The plots are shown in Figures 4 
to 9. To avoid awkward decimal values, units of the plots are deci- 
wavelength; hence 250 is 2.5 wavelengths, and 25 is a quarter of a 
wavelength. Recall that the beam center line is used as reference for 





Figure 5. Contours of Constant 
Phase Shift for 
0 = 18°. 



the phase shift. Reference to equation (51) shows that a positive 
phase shift is the case of more waves along the ray in question than 
along the ray coincident with the beam centerline. 

For the series of plots in Figures 4 to 9 the critical Mach 
number was selected. For a sphere this is M* = 0.587. The gap inte- 
gral is not included in the phase shift. The outer edge of each plot 
is for D/2R =0.5. Recall that D = beam diameter, and R = turret radius. 
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Figure 6. Contours of Constant 
Phase Shift for 
6 ’ = 36 ^ 




Figure 8. Contours of Constant 
Phase Shift for 
0' = 72°. 




Figure 7 . Contours of Constant 
Phase Shift for 
e' = 54°. 




Figure 9. Contours of Constant 
Phase Shift for 
6’ = 90°. 
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To provide insight to the magnitude of gap integral, Table II 



gives the maximum value for each of the plots. 



TABLE II. MAXIMUM VALUE OF GAP DTTEGRAL 



beam elevation angle maximum value of gap integral 

in deciwave lengths and location 
6 ' within beam 

(gap integral) a 



0° 


+ 36.8 




o 
00 
1 — I 


+ 65.4 


0 


36° 


+ 67.8 


0 


54° 


- 83.6 


180 


72° 


- 86.3 


120 


O 

o 

CD 


- 86.6 


90 



M = 0.586 for these calculations. Also D/2R = 0.5. 

oo 

* not a function of a. 



Evaluation of Phase Distortion as a Function of Distance Along the Beam 
An important question is the rate of decay of the integrand in 

equation ( 51 ) . Calculations were made of 

8(N - Nq) 

9(s/R) 

which is the phase distortion per unit distance along the beam. The 
phase distortion per unit distance is a function of distance along the 
beam. Results of the calculation are shown in Figure 10. The phase 
distortion is the area enclosed by one of the curves. Most of the phase 
distortion occurs within a length equal to one turret radius. 
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RATIO OF DISPLACEMENT ALONG BEAM 
TO TURRET RADIUS. Vp 

Figure 10 . Phase Distortion Per Unit ( s/R) as a Function 
of Distance Along the Beam for Rays at 
a = 0°, 60°, 120°, and 180° Within the Beam. 
Elevation Angle was 9 ' = 18° . 



Dependence of Phase Distortion on Flight Mach Number 

Another parameter of considerable interest is the dependence of 
phase distortion on Macli number. Qualitatively, isodensity contours at 
< M* have an appearance the same as for M*. Quantitatively, there 
are changes in both the shape of an isocontour of constant phase shift 
and the magnitude of the phase shift. Table III compares mgnitudes 
of distinctive features of the maps for two Mach numbers. Table III 
should be studied simultaneously with the appropriate figure for 9 ' . 

Phase distortion curves based solely on a first order calculation, 
i.e., using <('q in equation (9) only, should have the same shape for all 
M^. A calculation of phase distortion based on a second order solution, 
i.e., using both ((iq and in equation (9), should be a function of M^. 
Figure 11 shows the change in shape of a phase distortion curve with 
Mach number. The curves are drawn for G’ = 18°; M , either 0.300 or 
0.587; R'/R from 0 to 0.5 and for a = 0° and a = 180°. The phase dis- 
tortion has been normalized with the value at R’/R = 0.5 and a = 180°. 
Differences in the curves are due to Mach number effects. 
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TABLE III. COMPARISON OF MAGNITUDE OF PHASE SHIFT FOR TWO 



DIFEEI^ENT MACH NUMBERS 



0' 


feature of isocontour map 


value at M = 0.300 
00 


M = 0.1 
00 


0° 


maximum (|) at beam edge 


-41.5 


-160 


18° 


maximum <() at positive peak 


+15 


+ 60 




minimum ()> 


-69.6 


-290 


36° 


maximum ()) 


+44.2 


+240 




minimum (j) 


-80.3 


-345 


54° 


maximum <p 


+85.2 


+370 




minimum ()> 


-54.4 


-268 


72° 


maximum <|) 


+ 83.7 


+386 




minimum (j) at negative valley 


-14.1 


- 74 


90° 


maximum <() at a = 0° 


+49.1 


+ 223 




Figure 11. Phase Distortion for Two f-lach Numbers. 
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Interpretation of Results 



Phase shift for 6 ' = 0 as shown in Figure 4 consists of a series 
of concentric circles. Since the phase shift is negative, the beam 
is being focused. Concentric curves with negative phase distortion 
gives a positive focus. 

The contours of constant phase shift for 6 ' = 54° are almost 
straight lines. The algebraic signs of phase shift contours indicate 
the beam is being tilted in a direction opposite to the relative wind. 
The angular tilt is given rou^ly by 

(tilt angle) = |^:N-Nq)qo - ^ (52) 

Inserting the appropriate values (values for N-Ng were obtained from 
Table III) gives 



(tilt angle) 



(370 + 268) (3.8 x 10 ^m) 
(2) (.4572m) (100) 



26 .5 microradians 



Tilt of a beam does not cause a serious problem. 

These results for tilt angle can be compared with similar work by 

r V (12) 

Cook. 

( 19 ) 

Captain Cook used the following values in a sample calculation 
(notation has been adjusted to match this report) : 



R' = 0.25m A =10.6 micron 

oo 

R = 0.71m ><'/Poo =0-22 cm^/gm 

M=0.5 k'=2.7x 10"^ 

oo 

For 9' = 54°, a beam tilt of 63 microradians was calculated. When adjust 
ments are made for different parameters used to calculate the two results 
i.e., 63 and 26.5 microradians, the answers agree quite favorably. 

It is possible to determine an approximate focal length, F^. The 
subscript, a, indicates the focal length in a plane defined by a of 
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Figure 2. The appropriate equation is 



1 _ 

F ■ (AR*)2 
a 



(N-Nq) 



R- 



- 2(N-Nq) 



4 



(53) 



Equation (53) was evaluated for a = 0 and for the following values for 
other parameters 

= 1.23 X 10 ^ gm/cm^ 

X = 3.8 X 10"® m M = 0.587 

00 • 00 

AR'= 0.1524 m R = 0.9144 m 



(N-Nq) at R^ = 0.799 wavelengths 
(N-Nq) at Rj = 0.214 wavelengths 
0 ' = 0 ° 

A focal length of 16.5 km was found. 

Before proceeding, some additional information concerning equation 
(53) is appropriate. Equation (53) is a specialization of the equation 



1 

F ■ 3x2 

X 



where 



I = 



b 

/■ 



ds 



For the problem at hand equations (54) and (55) yield 



(54) 



(55) 



X 



K i f^2 ^0 ^ ^ f ^ ^ ^0 ^0 J j 

'^\j~ J J 7 



(56) 



where ~ p(^ 2 ^’ *^1 ~ p(x^),etc. It is obvious that the third integral 

( 12 ) 

IS zero. Equation (56) leads to equation (53). From Cook's report for 
the previously cited values for R, R' , etc. , a focal length of about 
6 km is calculated. Adjusting the value for values of R, R' , M , etc. , 
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used in the calculations of this report, good qualitative agreement 
is obtained. 

(13) 

Wolters and Laffay have calculated focal length for an 
aperture of a turret on a F-15. The turret geometry was a cylinder 
with axis parallel to the frees tream; radius of the cylinder was 0.5 m. 

When viewed from a position normal to the axis , a cylinder appears as a 
rectangle. The square corners were rounded with a radius of 0.3 m. 

Without rounding of the comers, the cylinder viewed from its end would 
appear as a disc with radius 0.5m. With the rounding of the comers , 
the disc is reduced to 0.2 m; the disc is normal to the freestream 
velocity vector. 

(13) 

Wolters and Laffay determined the focal length of the lens 
created by density variations within the external flow. The following 
parameters were used in their calculations: 

R' = 0.17 m X = 10.6 X 10“^ m 

00 

R=0.5m K*/p= 0.223 x 10~^ m^/kg 

c °° 

M = 0.587 0' = 0 

00 

altitude: 20,000 ft. 

For these conditions a focal length of 18.2 km was determined. Since 
the distortion varies linearly with the ambient density, a fact which 
follows from equation (42), the focal length can be adjusted to sea level 
by an adjustment to density. For the same conditions as above except at 
sea level, the focal length would be 8.2 km. A value 8.2 km is consistent 
with results obtained in this report, 16.5 km. 

Having compared in a rough way calculated results from references 
(12) and (13) with calculated results from this paper, the relative 
advantages and disadvantages of each model will now be compared. Reference 
(12) used 9I/9X and 9l/9y to predict tilt; 92l/9x^ [see equation (54)] and 
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32 l/ 9 y 2 were used to predict focal length. These were evaluated at the 
beam center. An analytical solution was obtained. Reference (12) also 
used a first order solution for the flow field, which is equation (16). 
This paper uses a second order solution as given by equation (41). Hence 
a more accurate flow field is used in this paper. Reference (12) yields 
a solution in simple form; this paper requires a desk top computer, eg. 
HP9830, or else a high degree of patience if a HP55 were to be used. 

In contrast to reference (13) which used a numerical solution to 
the flow field, the second order analytical solution gives the flow field 
to any degree of spatial resolution required. Numerical solutions give 
results in discrete steps. Based on Figure 5 of reference (13) the ratio 
of step size to cylinder radius is 0.22. Figure 10 of this paper shows 
that approximately one-half of the phase distortion is generated in a 
step size 0.22. Of course a smaller step size can be used. Numerical 
techniques have the advantages of providing solutions for complex turret 
geometry. This paper is applicable only to a hemispherical turret. 

Look at Figure 9 for which 6* = 90°. The contours have a symmetry 
about the a = 0°--180° line and about the a = 90°--270° line. The phase 
distortion has positive values which give a negative lens. The focal 
length for a = 0° is considerably less than for a = 90°. At 6' = 90°, 
the beam is defocused. 

It is apparent from Figures 4 to 9 that the phase distortion is a 
strong function of elevation angle 9 ' . Compensation for the phase dis- 
tortion is possible, at least partially. However, the shape of the pri- 
mary telescope mirror required for conpensation is a function of both 
elevation angle and Mach number. Such compensation should fall within 
the capability of adaptive optics techniques. 
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Both 



To discuss the altitude dependence, recall k = K'p^/pgj^. 

K* and Pgj^ are constants for a given wavelength and specified atmospheric 
conditions. Consequently, k varies directly with p^. Now examine equa- 
tion ( 51 ). Since k << 1 , the term 1 + < is essentially unity; it follows 
that N-Nq varies directly as k. Phase distortion changes as p^. Obviously 
as altitude increases, phase distortion decreases. 
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VIII. CONCLUDING REMARKS 



Compressible flow over a sphere has been obtained forming the basis 
for a model of a blunt hemispherical laser turret. Phase distortion has 
been calculated for several freestream Mach numbers. Phase distortion 
data are presented for M^ equal to the critical Mach number. Calculations 
are consistent with those of references (12) and (13). 

The phase distortion is a function of both elevation angle and Mach 
number. Hence optical compensation requires a variable shape primary 
mirror. 

The model of a hemispherical laser turret provides insight to the 
behavior of phase distortion as a function of I'fach number. See Figure 
11. Also refer to Table III. 

Another valuable aspect of the calculations is the decay of den- 
sity perturbation with distance along the beam. See Figure 10. It is 
apparent that the distortion is fairly localized. Within a distance of 
s/R = 1.0, over 99% of the distortion occurs. 

The formulation clearly shows the altitude dependence of phase 
distortion. 
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